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1. Introduction 

The development of the heavy quark or Isgur-Wise symmetry |IJ has stimulated a great 
deal of interest in studying [@H|3 properties of heavy baryons (i.e., those with the quark 
structure qqQ) in the bound state approach. In this picture the heavy baryon is treated 
as a heavy spin multiplet of mesons (with structure Qq) bound in the field of the nucleon 
(qqq) which itself emerges as a soliton configuration of light meson fields. This treatment 
is suggested by the l/N c expansion || of QCD. Recent reviews of the soliton approach to 
the light baryons are given in refs J9|— 112|| while the bound state treatment of the "light" 
hyperons is discussed in refs fl3], [T^| . 

A compelling feature of this approach is that it permits, in principle, an exact expansion 
of the heavy baryon properties in simultaneous powers of 1/M, 1/Nq and, since it is based 
on a chiral Lagrangian, number of derivatives acting on the light components of the heavy 
system. In practice there are obstacles related to the large number of unknown parameters 
which must be introduced. Rather than treating the light soliton in a model with many 
derivatives of the light pseudoscalar fields it turns out to be much more efficient to use the 
light vector mesons. Based on a model []15] of the light vector interactions with the heavy 
multiplet, the leading order (in the 1/Nq and 1/M expansions) heavy baryon mass splittings 
have been discussed |16|], obtaining satisfactory agreement with experiment. Actually the 
need for light vector mesons is not surprising since, in the soliton approach, they are necessary 
to explain, for example, the neutron-proton mass difference |T7]] and the nucleon axial singlet 
matrix element [18[ . 

In the present paper we focus our attention on the hyperfine splitting, which is of sublead- 
ing order both in 1/M and 1/Nq. This is a more complicated calculation and also involves 
using a cranking procedure |19] to obtain physical states which carry good spin and isospin 
quantum numbers. The first calculation of the heavy baryon hyperfine splitting in the per- 
turbative bound state framework was carried out by Jenkins and Manohar who got the 
formula 

mm) - m(s g) = < m < A) - rsM [ M - - Ml , (L1) 

\ Q) v QJ 4dF'(0) ' V ; 

where M* — M is the heavy vector-heavy pseudoscalar mass difference, d is the light 
pseudoscalar-heavy meson coupling constant and -^'(0) is the slope of the Skyrme "pro- 
file function" at the origin. This formula is obtained (see also section 5) by using the leading 
order in number of derivatives (zero) and leading order in 1/M heavy spin violation term. 
Therefore it is expected to provide the dominant contribution. Unfortunately, on evaluation, 
it is found to provide only a small portion of the experimental S*-S c masss difference. This 
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naturally suggests the need for including additional higher order in derivative heavy spin 
violation terms. However, there are many possible terms with unknown coefficients so that 
the systematic perturbative approach is not very predictive. 



To overcome this problem we employ a relativistic Lagrangian model |15j which uses 
ordinary heavy pseudoscalar and vector fields rather than the heavy "fluctuation" field mul- 
tiplet 0. This model reduces to the heavy multiplet approach in leading order and does 
not contain any new parameters. In a recent note |2(J we showed that such a model (con- 
sidered, for simplicity, to contain only light pseudoscalars; i.e., the light part is the original 
Skyrme model ||21|| ) yielded a "hidden" heavy spin violation which is not manifest from the 
form of the Lagrangian itself. This hidden part involves two derivatives and is actually more 
important numerically than the zero derivative "manifest" piece which leads to eq (|1 . 1| ) . 
However this new result is still not sufficient to bring the predicted S*-S c mass difference 
into agreement with experiment. The prediction for this difference is actually correlated to 
those for £ C -A C and A—N, the A - nucleon mass difference by [ |13| : 



m (E*) - m (E c ) = m (A) — m (N) — - [m (E c ) - m (A)] . (1.2) 

This formula depends only on the collective quantization procedure being used rather than 
the detailed structure of the model. If m (E c ) — m (A) and m (A) — m (N) are taken to agree 
with experiment, eq ( |1.2| ) predicts 41 MeV rather than the experimental value of 66 MeV. 
This means that it is impossible to exactly predict, in models of the present type, all three 
mass differences which appear in eq (|1.2|) . The goodness of the overall fit must be judged by 
comparing all three quantities with experiment. Our focus, of course, is the left hand side 
of eq ( |1.2| ) which is of order 1/M while the right hand side involves the difference of two 
order M° quantities. A similar calculation in the model with only light pseudoscalars was 



carried out by Oh and Park [ 22[ . However, they did not make a 1/M expansion in order to 
reveal the hidden violation terms. They also introduced a one-derivative "manifest" heavy 
spin violation term with a new relatively large unknown constant in order to improve the 
agreement with experiment. 

In the present paper we show that it is not necessary to introduce any new violation terms 
to agree with experiment if a chiral Lagrangian including light vectors is employed. Typical 
results are summarized, compared with experiment and compared with the Skyrme model for 
the light sector in Table |1.1| . A much more detailed discussion is given later in the text. We 
notice from the last row, that the model with light vectors gives a very satisfactory account 
of the E*-E c hyperfine splitting in contrast to the model without light vectors. There are 
also noticeable effects when the use of the heavy meson reduced mass is taken as a simple 
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Table 1.1: Typical results for the present model (including light vectors) com- 
pared with model with light pseudoscalars only ("Skyrme" column) and com- 
pared with experiment. No "manifest" heavy spin violation effects other than 
M* ^ M have been included. The column "present model + CM" simply takes 
into account recoil corrections by replacing the heavy meson mass by the reduced 



mass. A' c denotes a negative parity, spin 1/2 state. The quantity a in eqs Q2.6| ) 
and ( |2.7| ) was taken to be zero. All masses in MeV. 



mass difference 


expt. 


present model 


present model + CM 


Skyrme 


Ac-N 


1345 


1257 


1356 


1553 


A b - A c 


3356±50 


3164 


3285 


3215 


A' c - A c 


308 


249 


342 


208 


S c -A c 


168 


172 


158 


185 




66 


42 


63 
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approximation for kinematical corrections. Similarly, the first four rows of Table [TT] show 
that the other predictions of the model with light vectors agree well with experiment. Note 
that (see section 4) the predictions for mass differences are considered more reliable than 
those for the masses themselves. 

The present article is organized as follows. In section 2 we will review the classical, 
leading in 1/Nc, part of this calculation. This discussion includes both the light and heavy 
meson pieces of the Lagrangian. The emergence of bound state solutions is also explained 
in section 2. In section 3 we will describe the collective quantization in the framework 
of the cranking procedure for the bound states. Section 4 contains a detailed discussions 
of the numerical results. In section 5 we will discuss some new manifest contributions to 
the hyperfine splittings and the extension to different channels in the framework of the 
perturbation approach. We will conclude in section 6. The explicit expressions for the 
couplings between the bound heavy meson and the collective coordinates are listed in the 
appendices. 

2. The Model Lagrangian 

In this section we review the classical, i.e. leading order part in the 1/Nc expansion of 
the bound state description for the heavy baryons in the soliton picture. 

2.1. Light Mesons 
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For the sector of the model describing the light pseudoscalar and vector mesons we 
adopt the chirally invariant Lagrangian discussed in detail in the literature [p3j ^4]. This 
Lagrangian can be decomposed into a regular parity part 



tr 



U + [/t _ 2 ] _ tr [F^ (p) 0)] + rn^tr [R^] (2.1) 



and a part which contains the Levi-Civita tensor, e^ap. The action for the latter is most 
conveniently displayed with the help of differential forms p = p^dx^ 1 , etc. 



2N C 

157T 2 

+ (Tr 



Tr(p 5 
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t(Ti + ~l2)Rp A - ~-y 2 F(p)(pR - Rp) - 2ig\ l2 + 2 l3 )R i p 



3 w 2 ' 7 1 2 
In eqs (|2.1|) and ( |2.2|) we have introduced the abbreviations 



(2.2) 



(2.3) 



Here £ refers to a square root of the chiral field, i.e. U = £ 2 . Furthermore F^ v (p) = 
d^pv — d u pfj, — ig [p^pv] denotes the field tensor associated with the vector mesons p and 
u, which are combined in p^ = (u^l + p a ^r a ^j /2 when the reduction to two light flavors is 
made. The parameters g, 71, etc. can be determined (or at least constrained) from the study 
of decays of the light vector mesons such as p — ► 2n or uj — > 3n |24|] . 

The action for the light degrees of freedom (/ C s + r an ) contains static soliton solutions. 
The appropriate ansatze are 



cxp 



-r ■ rF(r 



uj\r} 



and Pi, a (r) 



G(r) 



e ija r j 



(2-4) 



9 gr 
while all other field components vanish. The resulting non-linear Euler-Lagrange equations 
for the radial functions F(r),u(r) and G(r) are solved numerically subject to the boundary 
conditions F(0) = —ir, uj'(0) = and G(0) = —2 while all fields vanish at radial infinity 
23. These boundary conditions are needed to obtain a consistent baryon number one 



configuration. 

2.2. The Relativistic Model for the Heavy Mesons 

In this subsection we present the relativistic Lagrangian, which describes the coupling 



between the light and heavy mesons [15] 



C 



H 



D^P (ZPp)t - t Q (Q^ - M 2 PP^ + M* 2 Q,Q 



2^ 



+2iMd (P P(t Qrt - Q^P ] 
2\f2icM f 



-e 

2 



my 



\2Q,F^ (p) Qt - -e^ 



Qv<xPnQ\ + QpVv {Q 



DpPF^ (p) Ql + Q a F,j, v (p) (D p Py 



(2.5) 
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Here we have allowed the mass M of the heavy pseudoscalar meson P to differ from the 
mass M* of the heavy vector meson Qy. Note that the heavy meson fields are conventionally 
defined as row vectors in isospin space. The covariant derivative introduces the additional 
parameter a: 

DyP ] = {dy.-iagpy.-i {I- a)vy)P ] = (dy - ivy - igaFLy) P ] , (2.6) 
DyQl = {dy-iVy-igaRy)Ql . (2.7) 

The covariant field tensor of the heavy vector meson is then defined as 

(Q,u) j = DyQl - D V Q\. (2.8) 



The coupling constants d, c and a, which appear in the Lagrangian ( |2.5| ), have still not been 
very accurately determined. In particular there is no direct experimental evidence for the 
value of a, which would be unity if a possible definition of light vector meson dominance 
for the electromagnetic form factors of the heavy mesons were to be adopted P5| . We will 
later adjust a to the spectrum of the heavy baryons. The other parameters in (|2.5|) will be 
taken |25| to be: 

d = 0.53 , c = 1.60 ; 
M = 1865MeV , M* = 2007 MeV , D - meson ; 
M = 5279MeV , M* = 5325 MeV, B - meson. (2.9) 

It should be stressed that the assumption of infinitely large masses for the heavy mesons 
has not been made in (|2.5| ). However, a model Lagrangian which was only required to exhibit 
the Lorentz and chiral invariances would be more general than the relativistic Lagrangian 
( p.5| ). The additional restrictions arise from the heavy quark transformation 

P' = e iMV - x P , Q' M = e iM * v - x Qy , (2.10) 

where the four-velocity V^ 1 characterizes the reference frame of the heavy quark. The heavy 
pseudoscalar and vector meson fields may then be combined in the heavy multiplet 

H=^{l + ~ Ul V»)(i lh P' + 1 v Q'u) and H = l0 H^ . (2.11) 
In the heavy quark limit, M = M* — > oo, the relativistic Lagrangian ( ^.5[ ) becomes 

jj£ H - iVTr {HDyR} - cfTr {Hj^S} - i^-Tr {H^ v F^(p)H} + ... . (2.12) 
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The ellipses indicate subleading pieces in 1/M. Actually the coefficients of the various 
Lorentz and chirally invariant pieces in the relativistic Lagrangian Q2.5|) have precisely been 
arranged to yield the spin-flavor symmetric model ( |2.12|) in the heavy quark limit [|15[| . 

2.3. Bound States 

Here we briefly review the origin of bound states in the S- and P-wave heavy meson chan- 
nels. These orbital angular momentum quantum numbers refer to those of the pseudoscalar 
component (P^) of the heavy meson multiplet (P\Ql). 

For the P-wave channel the appropriate ansatz reads 

Pt = -L$(r)f • rpe i£t , Q\ = -^ (r)pe lrf , (2.13) 

V47T V47T 



pe let . (2.14) 



'47T 

Note that here p refers to a properly normalized spinor which describes the isospin of the 
heavy meson multiplet. Similarly the ansatz for the S-wave is given by 

Pt = -L<S>( r )pe iet , Q\ = -L*„(r)f ■ rpe^ , (2.15) 

V47T V47T 



Q\ = -^[^ 1 (r)r i r-T + ^ 2 ( y r)rT-d i r}pe iet . (2.16) 

V47T 



It should be remarked that the isospin matrices, which multiply the isospinor p, have (since 
there are no unmatched indices) vanishing grand spin G, which is the vector sum of total 



spin and isospin. The above ansatze are substituted in the relativistic Lagrangian (2J5) and 
the resulting action functionals are listed in appendix A. The variation of these functionals 
yields the associated equations of motion. They are numerically integrated by adjusting 
the energy eigenvalue e so that continuous normalizable configurations are obtained*. This 
value of e directly yields the binding energy of the heavy mesons. The fact that we have 
U(r = 0) = — 1 at the spatial origin causes the angular barrier for the P-wave heavy meson 
to vanish while the S-wave acquires a finite one. As a result the P-wave heavy meson is 
more strongly bound. 

Finally we would like to mention the connection to the heavy quark limit. In that case 
the multiplet H is characterized by a single radial function |26[. This implies that in the 
limit M = M* — > oo the radial functions, which parametrize the bound heavy mesons, have 
to satisfy the linear relations 

^ = ^ 2 = -2$ (P-wave), (2.17) 



"The normalizability condition is quite restrictive. For example, it was shown in ref [|16[ to prohibit a 



''pentaquark" solution which would be extracted from eq (2.12) in the heavy quark limit. 
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ijr 1 = ^ 2 = $ (S - wave) , 



(2.18) 



together with \l/o — in both cases. Indeed the numerical solutions confirm these relations 
as the heavy meson masses approach infinity Hl6| . It should be remarked that commonly 
more than one bound state exists in each channel. Here we will concentrate mainly on the 
lowest one, which is characterized by the radial functions having no nodes away from the 
boundaries r = and r — > oo. However, in special instances we will also discuss the first 
radially excited state. 

2.4- Normalization 

As the relativistic Lagrangian (|2.5|) is bilinear in the heavy meson fields the resulting 
equations of motion are linear. Hence the overall magnitude of the solution is not fixed 
by the equation of motion. Nevertheless, the equations of motion for the heavy meson 
fields allow us to extract a metric for a scalar product between different bound states. In 
particular its diagonal elements serve to properly normalize the bound state wave-functions. 
The Lagrange function which results from substituting the ansdtze Q2.4Q and (|2.13|) - (|2.15|) 
may generally be written as 

L = -M d [F, G, u] + J e [F, G, u; $, %, * 2 ] p f (e)p(e) . (2.19) 



Here M c \ denotes the soliton mass |24]] whose minimum determines the light meson profiles 



F, G and to. The explicit expressions for the functionals I t are given in appendix A. The 
subscript refers to the explicit dependence on the energy eigenvalues. Upon canonical quan- 
tization the Fourier amplitudes p(e) and p'(e) are respectively elevated to annihilation and 
creation operators for a heavy meson bound state with the energy eigenvalue e. Demanding 
that each occupation of the bound state adds the amount |e| to the total energy yields the 
normalization condition]] 

[$,*o,*i, # 2 ] | = 1 (2.20) 

in addition to the canonical commutation relation [pi(e), pl(e')] = <%5 e ,e'- Note that for 
bound states the energy eigenvalues are discretized. For the P-wave channel we obtain the 



Hn the case that the explicit dependence on e is quadratic the proof is sketched in section 3 of ref 



A simple verification of eq (2.20) may be also obtained from its close connection to the Noether charge 
for the heavy quark number conservation. The latter is gotten by transforming the heavy fields as 
p _> e -<'iWp ) Q M ^ e-^^Qp and computing the quantity N = j d 3 x <5£(e~"' P, . . .) /6d ri | =flb =0 < 
However because the present ansatz is of the form P = P{x)e~ let we may equivalently compute this as 
N = J d 3 x SC I e -i (' 7+et )P (a;) ,...)/§€ .In the one heavy quark subspace this yields 1 = dl e /de. 

\ / rj—dor]—0 
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P - wave 

E = 1388 MeV 



S - wave 

£ = 1618 MeV 



















— % 




— - -p, 




— *2 



1.0 1.5 
r*m v 




Figure 2.1: The profile functions for the bound state wave-functions in the P- 
wave (left panel) and S-wave (right panel) channels. These functions are mea- 
sured in units of my = 773MeV. See text for the specification of the remaining 
parameters. 



normalization condition 



drr 2 {2 e 



-d 



a 



-sinF^i - -F'm 2 
r 2 



a 



2 

4\/2c 1 



-R a *o + ( e - f w ) 



[G (G + 2) *i - G'r^ 2 ] $ 



from eq (|2.20| 



gray v- 

For convenience we have employed the abbreviation R 



-- 1 (2.21) 
cosF - 1 + 



a (1 + G — cosF). Similarly for the S-wave channel the condition fl2.2U| ) yields 



drr 2 {2 e 



+4 



a 



a \ _ # a + 2 

-a; * 2 

2 J 2r 



*1 



# 2 - 2d 



4v/2c 1 



[G (G + 2) + 2G'r# 2 ] $ 



-sinF^i + 



1 . 



(2.22) 



gray r c 

The radial profiles associated with these normalizations are displayed in figure 2T1 for the 
choice a = —0.3. The parameters in the relativistic Lagrangian Q2.5Q have been set to the 
charm sector (|2.9|) . The parameters entering the light meson Lagrangian ( |2.1| , [2.2|) are given 
in eq JO]). 

3. Cranking the Bound Heavy Meson State 



It can easily be verified that the field configurations for both the light mesons 
and the heavy mesons (|2.13 )-( 2.16 ) are neither eigenfunctions of the spin- nor the isospin 
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generators. Hence these configurations do not possess the correct quantum numbers. In order 
to generate states which correspond to physical baryons a cranking procedure is employed. 
In the first step collective coordinates, which parametrize the (iso-) spin orientation of the 
meson configuration, are introduced. 

3.1. Collective Coordinates 

Time-dependent solutions to the equations of motion are required to obtain non-vanish- 
ing spin and isospin as the corresponding Noether charges. Unfortunately these solutions 
are unknown. Taking, however, into account that static rotations in coordinate and isospin 
spaces do not change the potential part of the Lagrange function, the assumption that these 
rotations are time-dependent seems to be a reasonable approximation. We therefore extend 
the soliton ansatz ( |2.4j ) by 

£ — A(*)£At(f) and Pfl — A(t) Plt A\t) . (3.1) 

Note that p M contains both isoscalar and isovector pieces. It should also be remarked that 
introducing only an isospin rotation as in eq ( |3.1| ) is sufficient because the hedgehog structure 
of the classical configuration ( |2.4| ) allows one to express a spatial rotation as one in isospin 
space. The time-dependence of the collective rotations is most conveniently parametrized 
by introducing angular velocities ft via 

A\t)jA{t)= l -r.a. (3.2) 

In the specific case that the Lagrangian contains terms which are linear in the time derivative 
as in eq (|2.2|) , additional field components are induced. For the light vector mesons these 
are linear in the angular velocities 

2 

LUi = -v?(r)e iiA; n i f fc and p k = fi(r)fi fc + £ 2 (r)r ■ fif fc . (3.3) 



Substituting the configurations (|3.1| ) and ( pT3|) into the light meson Lagrangian yields a term 



which is quadratic in the angular velocities. The constant of proportionality defines the 
moment of inertia a 2 [F, G, u; £ 1; £ 2 > V 9 ]- The induced radial functions <p(r), £i(r) and £z(r) 
are obtained from a variational approach to a 2 ||28||f[ The resulting equations of motion are 
coupled inhomogeneous differential equations with the classical fields F, G and u, which are 
fixed from extremizing the classical mass, acting as sources. Here it is worth mentioning 
that a 2 is of the order Nq- 



*The term a 2 refers to a frequently adopted notation for the moment of inertia and should not be confused 
with the coupling constant a in eq (127 
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Since the heavy mesons are also isospinors the collective rotation has to be applied as 
well. In analogy to eq (|3.1| ) we write 



P t_^ A(t)P t and Qt— ACOQj;. (3.4) 

Substituting the collectively rotating configurations into the total Lagrangian finally yields 

L P = -M cl + If )pt p + I a 2 ft 2 + l - Xp p^ ■ rp . (3.5) 

For convenience we have omitted the argument of the iso-spinor p. The classical mass M c \, 
which upon minimization provides the soliton profiles (|2.4 ), and the moment of inertia a 2 



are functionals of only the light meson fields. The quantity if 3 is given in eq and has 



already been employed to obtain the bound state P-wave profiles fl2Jj| ). The new quantity 



is the hyperfine parameter xp whose explicit expression is displayed in appendix B (|B.1|) . 



3.2. Quantization of the Collective Coordinates 

Here we will discuss how the canonical quantization of the collective coordinates A leads 
to states which may be identified with physical baryons. In order to construct Noether 
charges we first have to consider the variation of the fields under infinitesimal symmetry 
transformations. For the isospin transformation we observe 



.Ti 

, i— 
2 



+ (3.6) 

Here refers to any of the iso-rotating meson fields and the ellipses represent terms, which 
are subleading in 1/Nc, as e.g. time derivatives of the angular velocities which might arise 
from eq (|3.3j ). Furthermore Dij(A) = (1/2) trfjiArjA') denotes the adjoint representation 
of the collective rotations A. From eq ( |3.6|) we conclude that the total isospin is related to 
the derivative of the Lagrange function with respect to the angular velocities 

dL 

Next we note that the total spin operator J enters the grand spin operator G in the labo- 
ratory frame via 

Gi = Ji + D^ 1 (A)Ij — Ji— Jf . (3.8) 

For convenience we have defined the spin carried by the soliton «J so1 = dLp/dVt. As a 
consequence of the relation (|3.7|) , its absolute value is identical to that of the isospin, i.e. 
(«7 so1 ) 2 = I 2 = 1(1 + 1). By construction the light meson fields do not contribute to G. Even 
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more importantly and as has been noted before, the pieces of the heavy meson wave-functions 
( p.4|) , which multiply the spinor Ap, have zero grand spin too. Using the normalization 
condition (|2.21|) one therefore ends up with 

G = -p^p . (3.9) 

This relation will be helpful because it relates the operator multiplying the hyperfme param- 
eter in the collective Lagrangian ( |3.5|) to the spin and isospin operators. The collective piece 
of the Hamiltonian is obtained from the Legendre transformation 

1 



H%° 11 = ft ■ J so1 - L coU 



J so1 + XpG~ , (3.10) 



Here Lp 11 refers to the fl dependent terms in eq ( |3.5| ). Finally the mass formula for an even 
parity baryon with a single heavy quark becomes 

M P = Me + \e P \ + ^ [ XP J(J + 1) + (1 - XpW + 1)] , (3-11) 

where contributions of 0(xp), which apparently are quartic in the heavy meson wave- 
function, have been omitted for consistency because terms of that order have been excluded 
from the very beginning. Also the omitted terms are sub leading in the 1/M expansion 



since xp goes as 1/M, cf. figure fO] and Appendix B. In addition, the operator contained 
in the omitted term, (p^rpj ■ (pWp), does not contribute to the hyperfme splitting. The 
reason is that from canonical commutation relations for the components of the isospinor p, 
Pi, plj = 5ij, this operator is shown to be Nq (Nq + 2), where Nq is the occupation number 
for the heavy meson bound state [13|. Hence this term contains neither the spin nor the 
isospin quantum numbers. 



From eq ( |3.11|) we recognize that the spin degeneracy between baryons containing a heavy 
quark vanishes in the heavy quark limit because xp approaches zero. Of course, this result 
is a direct consequence of the spin-flavor symmetry and would not have come out in case 
the various Lorentz and chirally invariant terms in eq (|2.5|) had been chosen arbitrarily. 



3.3 The Odd Parity State 

The S- and P-wave heavy channels decouple because they have opposite parity. Therefore 
the quantization of the S-wave bound state may be considered independently from the P- 
wave case, which has been discussed in the preceding section. The calculation, which proceeds 
along the lines of the one discussed in subsection 3.1, yields 

L = -M cl + If Vp + ^a 2 n 2 + ^x 5 p f O • rp . (3.12) 
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Here, of course, the spinor p corresponds to the one of the bound heavy meson in the S-wave 
channel. The explicit expression for the corresponding hyperfine parameter xs is given in 
eq ( B.3| ) in Appendix B. 



We may apply the same quantization procedure as for the P-wave. This results in the 
mass formula 

M s = M cl + \e s \ + ^2 iXsJ(J + 1) + (1 - XsW + 1)] (3.13) 

for baryons constructed as a light baryon and a single occupation of the bound state for the 
heavy meson being in the S-wave channel. 

Let us add a brief comment on the 1/Nc dependences in eqs an d ( P-13|) . The 

classical mass is 0{Nc) while in leading order the binding energies are 0(N ). As already 
noted above, the moment of inertia is O(Nc) while x ~ ®{N C ). Therefore the hyperfme 
splitting is not only subleading in the heavy quark limit but also in the 1/Nc expansion. 

4- Numerical Results 

In this section we will discuss the numerical results obtained for the masses of the heavy 
baryons in the model discussed above. In particular we will concentrate on the spin and 
isospin splitting in the realistic case of finite heavy meson masses Q2.9| ). It should be noted 
that sizable quantum corrections occur for the classical soliton mass M c \ [p9fl . It seems that 
these corrections are (approximately) equal for all baryons. Hence we will only consider 
mass differences between various baryons. In that case the absolute value of the classical 
mass M c \ is redundant. The parameters in the light sector cannot completely be determined 
from properties of the corresponding mesons. The remaining (limited) parameter space is, 
however, more than fully constrained by a best fit to the mass differences of the low-lying 
| + and | + baryons in the light sector. This yields: 

# = 5.57, my = 773 MeV 

7i = 0.3, 72 = 1.8, 73 = 1-2. (4.1) 



The resulting mass differences for the light baryons all agree within about 10% PC||. The 



corresponding moment of inertia is a 2 = 5.00GeV . In eqs ( 3.11 ) and ( 3.13 ) 1/a 2 enters as 



the coefficient of those terms which determine the hyperfine splitting. Hence a fine-tuning 
of the parameters ( |4.1| ) to e.g. a 2 = 5.11GeV _1 , which exactly reproduces the A-nucleon 
mass difference, has only negligible effects on the predicted hyperfine splittings. 

Before discussing the implications for the physical parameter results ( [2.9|) we would like 
to comment on the heavy limit behavior of the hyperfine splitting parameters xp an d Xs- F° r 



13 



Hyperfine Splitting 




15.0 20.0 25.0 30.0 

M=M* (GeV) 

Figure 4.1: The hyperfine splitting parameters xp an d Xs a functions of identical 
heavy meson masses for different values of the coupling constant a. 



this purpose we have plotted these quantities as functions of M — M* in figure |4.1| . We see 
that for both channels the splitting parameters decrease when the heavy limit is approached. 
In the appendix we show that the leading order term in the heavy quark expansion indeed 
is proportional to 1/M. Clearly the hyperfine parameter in the S-wave channel decreases 
somewhat more quickly with the heavy meson mass than in the P-wave channel. 

For fixed M the hyperfine parameters in the two channels behave oppositely with regard 
to the undetermined coupling constant a: While xp decreases when a becomes larger, 
Xs increases. In ref |20 we have shown that a major fraction of the P-wave hyperfine 
constant is due to terms in the relativistic Lagrangian ( |2.5| ) which do not manifestly break 
the heavy spin symmetry rather than to terms, which explicitly break this symmetry; as 
for example M ^ M* . For a quantitative discussion of this hidden contribution we perform 
the calculation using identical masses from the charm sector i.e. M = M* = 1.865GeV and 
furthermore a = 0.3. We take all other parameters as in eq (|2.9|). This results in xp — 0.080. 
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From table |4.1| we recognize that this is about 80% of the value obtained using the physical 
masses M = 1.865GeV and M* = 2.007GeV. In the case of the S-wave the hidden piece is 
even more dominant. For the symmetric choice of the mass parameters one finds xs — 0.175 



which is more than 90% of the value displayed in table |4.1| for a = 0.3. It is also interesting 



to compare figure 4.1 with the corresponding curve in figure 1 of ref |^0fl , pertaining to the 
model without light vectors. This makes it clear that the light vector model predicts a 
substantially larger xp- 

In section 5, some more discussion of the "hidden" hyperfine splitting terms is given. In 
addition, three more "manifest" heavy spin symmetry violating terms associated with the 
relativistic Lagrangian ( [2.5|) are treated in the perturbative expansion. Since the manifest 
(M* — M) contribution is relatively small it is reasonable to expect that the others will be 
small too. 

Let us next discuss the spectrum of the baryons containing a single heavy quark. For 
this case we assume the realistic masses as in eq ( |2.9|) . In table |4.1| the numerical results for 
the lowest S- and P-wave bound states in the charm sector are displayed. As already noted 



in ref |16| the binding energy 

u = M-\e\ (4.2) 

decreases with growing coupling constant a. This is the case for both the P and S-wave 
channels. For M — > oo the heavy limit |26] 

u — 3 -dF>(0) + ^G"(0) - ^(0) (4.3) 
2 gmy 2 



will be attained^. As in the discussion of figure |]T we see that the hyperfine parameters in 



these two channels behave oppositely as functions of a. Here we have chosen to measure the 
mass differences with respect to the lightest charmed baryon, A c . Hence the mass differences 
with respect to E c and S* directly reflect the a-dependence of hyperfine parameter xp while 
the corresponding dependence of the binding energy ujp can be extracted from the splitting 
relative to the nucleon. In addition the splitting with respect to the negative parity charmed 
baryons reflects the a-dependence of the S-wave channel binding energy u>s- Finally the 
mass difference to At contains the energy eigenvalues and hyperfine parameters computed 
with the B and B* meson masses in eq ( |2.9|) . 

While the mass difference to the nucleon is improved with a positive value for a, the 
agreement for the mass differences between the heavy baryons slightly deteriorates when in- 



Note that the conventions in ref 1 26 differ from the present ones as explained in Appendix B of ref 
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Table 4.1: Parameters for heavy baryons and mass differences with respect to A c . 
Primes indicate negative parity baryons, i. e. S-wave bound states. All energies 
are in MeV. 
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0.346 


s c 
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415 


413 


408 


402 
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? 


335 


K* 


468 


467 


464 


461 


458 


? 


437 


N 


-1237 


-1257 


-1278 


-1299 


-1321 


-1345 


-1553 


A b 


3160 


3164 


3167 


3170 


3173 


3356 ±50 


3215 



creasing this parameter. Nevertheless, fair agreement with the experimental data is achieved 
for quite a range of a. 

Table |4.1| also contains the model predictions when the background soliton is taken from 
the basic Skyrme model pi , |19|| which does not include the light vector mesons. Here we 
have adjusted the only free parameter (eskyrme = 4.25) to reproduce the A-nucleon mass 
difference. From the A c -nucleon mass difference we observe that in comparison with the 
nucleon the masses of the heavy baryons are predicted about 200MeV too large. This 
confirms the above statement that the spectra of both the light and the heavy baryons can 
only be reasonably reproduced when light vector mesons are included. This conclusion can 



already be drawn from the too small binding energies [16|. The hyperfine corrections make 
only minor changes in the A;,-A c splitting. 



In table |4.2| we display the analogous predictions for the bottom sector. According to 
the heavy spin symmetry the binding energies of the P- and S-wave channels approach each 
other. Hence the mass differences between the even and odd parity baryons containing a 



bottom quark correspondingly decrease. As was already inferred from figure [0] we confirm 
upon comparison with table f|4] that xp decreases less quickly with the heavy meson mass 
than xs- Except for A& no empirical data for the masses of these baryons are known at 
present. These results for the mass differences among the bottom baryons are predictions 
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Table 4.2: Parameters for heavy baryons and mass differences with respect to 
A;,. Primes indicate negative parity baryons, i.e. S-wave bound states. All 
energies are in MeV. The empirical value for the relative position of the nucleon 
is 4701±50MeV [31]. 
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Yp 
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u s 
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Xs 
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189 
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K 
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164 
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363 


359 


358 
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^b 


375 


373 


371 


369 


367 


N 


-4397 


-4422 


-4446 


-4471 


-4494 



of the model which can, in the future, be compared with experiment. As could have been 
inferred from the next to last row in table |4.1| the absolute position of the bottom multiplet 
is about 200 ± 50MeV too low. On the absolute scale this apparently is only a 5% deviation 
from the data. Certainly a larger value a m 1, which corresponds to a model for light vector 
resonance dominance of the heavy meson form factor |25[] , would yield an excellent agreement 
for the mass difference between A;, and the nucleon. On the other hand such a choice would 
slightly spoil the nice picture for the charm multiplet. 

In table |4.3| we list the numerical results for baryons constructed from the first radially 
excited P-wave bound state in the charm sector. The particle data group (PDG) lists 



an excited A c (2625), about 340MeV above the A c [31|, although this state is more likely, 
according to the PDG, to have J p = | . In the bound state picture this would require 
a D-wave (G = 3/2) heavy meson bound to the soliton. Here we have not discussed that 
channel. A preliminary discussion in the perturbative approach is given in the next section. 

The preceding calculations are based on the Nq — > oo limit in which the nucleon is 
infinitely heavy. From a common sense point of view this is peculiar since the nucleon is 
actually lighter than the heavy mesons being bound to it in the model. Hence, for comparison 
with experiment it is desirable to estimate kinematic effects associated with the nucleon's 
motion. These are expected |26| to lower the binding energy of the heavy baryons which 
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Table 4.3: Parameters for radially excited heavy baryons and mass differences 
with respect to A c . A tilde refers to a radially excited P-wave baryon. All 
energies are in MeV. 



a 


-0.1 


0.0 


0.1 


0.2 


0.3 


UJp 


125 


113 


101 


90 


79 


Xp 


0.020 


0.014 


0.007 


0.000 


-0.007 


A c 


429 


419 


409 


398 


388 




625 


618 


607 


598 


589 



have up to now come out too high (see A c -N mass difference in Table 4.1, for example.). In 
order to estimate these kinematical effects in the bound state approach we have substituted 
the reduced masses 

1 1 

M 



M, 



cl 



— and 

M M* 



1 1 



(4.4) 



into the bound state equations. In a non-relativistic treatment this corresponds to the 



elimination of the center of mass motion | 26fl . The results for the spectrum of the heavy 
baryons obtained from the replacement §L 



are in displayed in table |47|. Again we consider 
a as a free parameter. We notice that there is a remarkable improvement in the prediction for 
the A;, mass, which was previously the worst one. The changes in some of the mass parameters 
can approximately be compensated by a suitable re-adjustment of a. For a ~ 0.0 to —0.4 the 
agreement with the existing data is quite reasonable. When using the reduced meson masses 
the S c baryon is always predicted a bit too light while it is too heavy when the physical 
masses are substituted in the bound state equation. For A^, the situation is opposite. While 
the use of the physical meson masses gives too small a mass, the substitution of the reduced 
masses gives too large a prediction for the mass of this baryon. These results indicate that 
kinematical corrections are indeed important. It should, however, be remarked that with this 
replacement the heavy quark limit \ —>■ cannot be attained because the mass parameters in 
the bound state equations will remain finite and hence the linear relations Q2.17 ) and ( |2.18| ) 
will not be satisfied. However, these relations are essential to verify the limit x ~~ * for 
M = M* — > oo. Nevertheless we think that the replacement (4^) provides sensible insight 
in the relevance of kinematical corrections. 

It is interesting to see how far the heavy quark approach can be pushed to lighter quarks. 
To answer this question we have considered the strange quark. In the corresponding kaon 
sector the P-wave is only very loosely bound when the physical masses are substituted. On 



the other hand sizable binding energies are obtained when the reduced masses are used [16 
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Table 4.4: Parameters for heavy baryons and mass differences with respect to A c . 
Primes indicate negative parity states, i.e. S-wave bound states. All energies 
are in MeV. In this calculation the reduced masses (4.4) enter the bound state 
equations from which the binding energies are extracted. The physical meson 
masses 1865MeV and 5279MeV are used when computing the mass differences 
to the nucleon and the A& from these binding energies. Radially excited states 
are omitted because they are only very loosely bound, if at all. The empirical 
data are taken from the PDG [31], see also [32]. 
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Table 4.5: Same as table 4.4 for even parity baryons in the kaon sector. 
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This behavior is somewhat different from the charm and bottom sector and can be understood 



by noting that the difference M* — M is considerably reduced when using (|4.4| ). In the heavy 
sectors (charm and bottom) this difference is small in any event. The resulting spectrum for 
the strange baryons is shown in table ^5|. The comparison with the experimental data shows 
that even the use of the reduced masses does not provide sufficient binding. In the S-wave 
channel the situation is worse, even when the reduced masses are substituted bound states 
are not detected unless a < —1.0. The failure of the present approach in the strange sector 
strongly suggests that for these baryons a chirally invariant set-up [|12|] is more appropriate. 

5. The Perturbative Approach 

The perturbative approach can illuminate several aspects of the hyperfine splitting prob- 
lem. This is due to the heavy quark symmetry which is naturally exploited by making an 
expansion in powers of 1/M using the heavy field formalism. Our starting Lagrangian ( |2.5| ) 
has been set up in such a way as to yield a heavy quark symmetric result as M —>■ oo 
when M = M* is assumed, cf. eq (|2.12|) . The perturbative 1/M expansion is more general 
(presumably exact) but less predictive. Thus the 1/M expansion provides a useful calibra- 
tion in the large M limit. Since it deals with perturbation matrix elements it provides us 
with a convenient classification of the various sources of hyperfine splitting. The method is 
also advantageous in that it can be extended, without too much algebraic work, to different 
channels of interest. On the other hand, once the particular channels of interest are settled 
on, it is clearly more convenient to employ the exact numerical solution, which efficiently 
sums up a class of 1/M corrections. 

The leading order Lagrangian ( |2.12j ) can be supplemented by terms which manifestly 
break the heavy quark symmetry to leading order (M° with the present normalization) as 
follows: 



M 



-C 



M-M*. 



H 



-Tr 



+ 



{d - d') 



Tr 
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,/2(r — d) r i r n 

-i— -Tr ltllv HF^{p)H + aV p Tr Ha^ (gpp - vp) Ha^ (5.1) 

my J 

Here the (M — M*) term measures the heavy spin violation due to the heavy pseudoscalar - 
heavy vector mass difference. The (d — d') term measures the heavy spin violation induced 
by choosing different coefficients for the fifth and sixth terms in eq (|2.5| ) (or see eq (2) in 
ref ||20|| ) while the (c — d) term corresponds to choosing different coefficients for the last and 
next-to last terms in eq (|2.5| ). Finally the at term corresponds to the leading term obtained 
by using different values of a in eqs (fO| and fl2~77|) . Note that (M — M*), (d - df), (c - d) 
and a all behave as 1/M. 

In addition to the terms in eq ( |5 . 1| ) , which manifestly break the heavy quark symmetry, 
there are, in fact, "hidden" violation terms contained in eq ( |2.5|) . The explicit expression for 
the hidden terms in the model without light vectors is given in eq (11) of ref [p0| . These were 



shown to exist (for the model without light vectors) in ref [20[] and arise from performing a 
detailed 1/M expansion of the relativistic Lagrangian. In the sense of the chiral expansion 
these terms carry two derivatives, but nevertheless turn out to be very important numerically 
for the case considered. In the previous section the numerical study has confirmed that this 
is also true when light vector mesons are included. It was shown (cf. fig 2 of ref ]20|]) that 
the dependence on d of the hyperfine splitting computed from these hidden terms using the 
perturbative approach generally matched the exact numerical calculation. Hence we shall 
not explicitly isolate the extra hidden terms due to the addition of the light vectors but shall 
content ourselves with the numerical treatment given in the preceding section. 

In the perturbative approach the collective Lagrangian involving the variable A(t) is 
obtained by substituting 

H(x,t) = A(t)7T c (x) , (5.2) 

where H c (x) is the heavy meson bound-state wave function, into the heavy field Lagrangian. 
Clearly this is the analog of the replacement ( |3.4| ). (The treatment of the chiral Lagrangian 
of the light pseudoscalars and vectors is the same as in section 3.) The bound-state wave- 
function is conveniently presented in the rest frame, = (1, 0), where 

*-(;:)• 

with a, I, h representing respectively the isospin, light spin and heavy spin bivalent indices. 
Due to the hedgehog structure of the soliton profiles, the calculation is simplified if we deal 
with a reduced wave-function obtained after removing the factor r ■ r: 

hl = %(r.r) ad ^ hl (5.4) 
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whereQ u(r) is a radial wave function, assumed to be very sharply peaked near r = for 
large M. In the leading order of 1/M there is no violation of the heavy quark symmetry and 
we may perform a partial wave analysis of ipdi,h 



ipdi,h(9,93,r,k) = J2 C r£k 3 -,g 3 Y r 3 tdi(k } k 3 )xh ■ 



(5.5) 



Here Y^ 3 stands for the standard spherical harmonics representing orbital angular momentum 
r and C denotes ordinary Clebsch-Gordan coefficients. The heavy spinor Xh is trivially 
factored in this expression as a manifestation of the heavy quark symmetry. Furthermore 
£,di{k, k 3 ) represents a wave-function in which the light spin and isospin are added vectorially 
to give K = /nght + Slight with eigenvalues K 2 = k(k + 1). The total "light grand spin" 



g = r + K 



(5.6) 



is the significant quantity in the heavy limit. The dynamics of the model dictates that the 
bound-states occur for k = 0, in which case £^(0,0) = e^/\/2. The bound-state wave- 
function simply is 



^(0,0,0,0) 



1 



- e dlXh 



57T 



The k = 1 unbound wave-function with no orbital excitation (r = 0) is 

1 



^dlO;93)Xh ■ 



(5.7) 



(5. 



When violations of the heavy quark symmetry are included, g is no longer a good quantum 
number. We define the grand spin, which is a good quantum number, as, 



G' — g + Shenvy ■ 

In the notation of eqs ( |5.7| ) and (|5.8|) we have the grand spin eigenstates 



^\(G' = G' 3 = 1/2) 
V&G" = G' 3 = 1/2) 



1 



--£dA 



2h 



57T 



(5.9) 

(5.10) 
(5.11) 



Note that in eq ( |5.10| ) the G ; 3 = +1/2 wave function is #2h since the index 2 corresponds to 
+ 1/2 for the anti-quark wave-function. The two states ( |5.10 ) and ( |5.11 ) differ with respect 
to their g and K labels. 



*We have removed a factor of 1/V47T compared to ref |2C(], since it is now carried by the spherical harmonic 
in eq (|5.5|). 
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Now let us consider the potential for the bound-state wave-function in the presence of 
the first heavy quark symmetry violating term in eq (|5.1| ). Substituting the G"-eigenstates 
^M 1 ) and ip^ from eqs ( |5.10 ) and ( |5.11D into eq (|2.12|) and the first term of eq (|5.1| ) yields, 
after a spatial integration, the potential matrix in the ip^-ip^ space: 

v = _d_m(3 \ + M-M 1 ( » V5 




where F(r) is defined in eq ( [2.4|) and F' = dF/dr. The first matrix shows that ijj^ is bound 
while ip^ is unbound in the heavy spin limit. Since the second matrix gives mixing between 
ipW and the latter must be included in the presence of effects which break the heavy 
quark symmetry. The diagonalized bound wave function is seen to be 

m V3M-M* m 

** ' - -^fw 4 ■ (5 ' 13) 

This is the proper wave-function to be "cranked" in order to generate the heavy spin viola- 
tion. Using it in eq ( |5.2j ), which is then substituted into the a = limit of the first term of 
eq ( [2.12 ), contributes a term in the collective Lagrangian 

v M*-M . . 

2^3 where X = j^- (5.14) 

By using the Wigner-Eckart theorem we may express this for states of either G' 3 as the 
matrix element of the operator • G' . For convenience we have chosen to consider our 
wave-function as representing the conjugate particle in eq (|5.4|) . Hence the matrix element 
of G' in this section differs by a minus sign from that of G defined in eq ( |3.8|) . The latter 
is the appropriate one when we form the total heavy baryon spin J = G + J so \ with 
Jf 31 = (dL/dfli). Then the collective Lagrangian, L co u may be written (see section 3) 

L coll = ^a 2 Q 2 - X £l G (5.15) 

which again leads to the Hamiltonian (|3.10|) . Substituting a 2 = (3/2) [m(A) — m(N)} in eq 
( p.lO| ) we get the well known formula, cf. eq ([1 .If) 

m(Z* Q ) - m(E Q ) = [m(A) - m(N)] X • (5.16) 

The purpose in deriving this again was to explain the perturbative method and our notation. 

Next we shall give some new perturbative "manifest" contributions to \ from eq (|5.1| ). 
When all these terms are included the potential V in eq ( 5.12|) is modified so that the properly 
diagonalized wave-function which replaces eq (|5.13|) becomes 



+ e^ (2) , (5.17) 



23 



with 

_^ ( M _ M *) + f (d - d') ng) + ^3^(0) - VQ C -^^ 

2dF , {0) + 2_V2cGm 

y ' grriv 



There are two types of contribution to \. The first type is analogous to eq (|5.14j ) and arises 



when eq ( |5.17|) is cranked and substituted into eq ( |2. 12| ) . The second type is obtained by 



substituting the leading order wave function ipW into the (c — d) and a terms in eq ( |5.1| ) 
The complete expression for \ resulting from the "manifest" heavy spin violation is 



X 



+ ? [8 -20(6(0)- 6(0))] -W2— V "(0) . (5.19) 

3 my 



The quantities £i(0), £2(0) and ^"(0) are defined in eq (|3.3| ). This formula may be useful for 
quickly estimating the effects of heavy spin violation in the coupling constants, which were 
not explicitly given in the previous discussion. Unfortunately there is no determination of 
the magnitude of these effects from the mesonic sector at present. In ref |5D[ the discussion 
of the "hidden" heavy contributions to \ was given for the Lagrangian with only light 
pseudoscalars. 

The hyperfine splitting just discussed is for the ground state or P-wave heavy baryons. 
It is of some interest to briefly consider the negative parity heavy baryons with one unit of 
orbital excitation. In the heavy spin limit these bound states correspond to the r = 1 and 
k = choice in eq (|5.5| ): 

V^(l,£3,l,0) = ^|yfxh. (5.20) 

The spin, Jii g ht of the "light cloud" part of the heavy baryon is gotten by adding this g = 1 
piece to the soliton spin «J so1 . For the 1 = (which implies J so1 = 0) heavy baryons one finds 
flight = 1 and the degenerate multiplet 

{A' Q (l/2),A' Q (3/2)} . (5.21) 

For the / = J so1 = 1 heavy baryons, Jn g ht can be either 0, 1 or 2 and we find the degenerate 
heavy multiplets 

Sq(1/2) , 
{£' Q (l/2),£' Q (3/2)} , 

{£' Q (3/2) , £' Q (5/2) } . (5.22) 
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In general, the situation is even more complicated and further discussion will be given else- 
where. At present there are experimental candidates |3l| for a negative parity spin 1/2 baryon 
A' c at 2593.6 ± 1.0 MeV and a negative parity spin 3/2 baryon A' c at 2626.4 ± 0.9 MeV. 

Since experimental information is available, it is especially interesting to consider the 
splitting between the two A'q states in eq ( 5.21| ). This splitting stems from the violation of 



the heavy quark symmetry. For the Aq type states the total spin coincides with the grand 



spin G so that eq (|5.21|) may be alternatively considered a G = 1/2, G = 3/2 multiplet. 
Since the good quantum number is G, we may in general expect the hyperfme parameter x 
to depend on G. The collective Hamiltonian takes the form 

rsol - x 



J sol + XG G i 

tfcoii = ^ • (5.23) 

On general grounds we see that for the case of the Aq's the collective Hamiltonian contribu- 
tion to the hyperfme splitting will be suppressed. Setting J so1 = in eq ( |5.23| ) shows that the 



hyperfme splitting is of order (% 2 ) or equivalently of order (1/M 2 ). Unlike the ground state 
which involves only the G = 1/2 P-wave channel, there is another possibility for hyperfme 
splitting here. It is allowed for the G = 1/2 and G = 3/2 bound state energies to differ from 
each other. In the Lagrangian with only light pseudoscalars this does not happen and the 
Aq(1/2) — Aq(3/2) splitting is of order 1/M 2 . However when light vectors are added, there 
are "hidden" 1/M terms, which violate the heavy quark symmetry as e.g. 



i Tr 



a a ^H lu F^{p)D a H + h.c. . (5.24) 



This term is likely to generate splitting for the multiplet (|5.21|) to order 1/M by giving 
different binding energies to the G = 1/2 and G = 3/2 channels. It would be very interesting 
to investigate this in more detail. 

Finally, we add a remark concerning an amusing conceptual feature in the computation 
of hyperfme splitting among the five Sq's in eq fl5.22|) . The total angular momentum of each 
state is given by 

G 

J = £^£g +S he , V y , (5.25) 

J light 

where we are now considering each operator to be acting on the wave-function rather than its 
complex conjugate. We have illustrated two different intermediate angular momenta which 
can alternatively be used to label the final state. If «/n g ht is used, we get the heavy-spin 
multiplets in eq (|5.22|) . On the other hand, when the hyperfme splitting is turned on, the 



choice G is convenient, because it remains a good quantum number. According to the laws 
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of quantum mechanics, we cannot simultaneously use both to specify the states, since the 
commutator 

flight , G 2 = 4iJlighf (SheavyXtf) (5.26) 

is generally non-vanishing. This means that we cannot uniquely trace the splitting of, say, 
the {s / Q (l/2),S' Q (3/2)} heavy multiplet in eq ( |5.22|) , as hyperfine splitting interactions are 
turned on. Physically, this causes a mixing between the Sq's of the same spin. Rather, we 
must look at the whole pattern of the five masses. On the other hand, the problem simplifies 
for the computation of the ground state hyperfine splitting in eq ( |5.16| ). In that case the 
bound state wave function is characterized by g = 0. Thus the commutator in eq (|5.26|) 
vanishes, and it is "trivially" possible to track the hyperfine splitting as a mass difference. 

6. Discussion and Conclusions 

In the framework of the bound state approach we have studied the hyperfine splitting 
for baryons containing a heavy quark. In this approach a heavy baryon is constructed from 
a heavy meson configuration bound in the background of a (chiral) soliton. Here we have 
limited ourselves to heavy mesons in the S- and P wave channels, which exhibit the strongest 
binding. The study has been motivated by the earlier observation that light vector meson 
fields are required in the soliton configuration in order to reasonably describe the spectra of 
both the light and the heavy baryons when all available information on coupling constants 
of the elementary mesons is incorporated. The inclusion of light vector mesons causes some 
technical difficulties because field components which vanish classically are induced when 
time-dependent collective coordinates are introduced in order to generate states with good 
spin and isospin from the soliton. One might argue that the better agreement in the vector 
meson model is due to an additional parameter a; however, we have observed that the 
agreement is achieved for quite a wide range of this parameter. In fact the binding energies 
vary by only about lOOMeV in the range —0.1 < a < 0.3. On the other hand the discrepancy 
between the empirical data and the predictions obtained from the Skyrme model soliton is 
about twice as large. In addition the vector meson model reasonably reproduces the relative 
(to the nucleon) masses for both the charm and the bottom sector. Furthermore the mass 
difference within a given heavy multiplet, i.e. the hyperfine splitting, has turned out not to 
be very sensitive to that parameter either. 

We also have estimated kinematical corrections by substituting the reduced masses. The 
comparison with the empirical data has certainly indicated that these corrections are im- 
portant. This simple non-relativistic substitution fails, however, to satisfy the heavy quark 
limit result, which states that the hyperfine splitting should vanish for infinitely large quark 
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masses. It thus seems interesting to further explore the kinematical corrections. 



As an extension of earlier work |20] we have illuminated the systematics of the 1/M 
expansion of the hyperfme splitting. The main conclusion is that the dominant contribution 
stems from terms in the relativistic Lagrangian which do not manifestly break the heavy 
quark symmetry. 

We have furthermore observed that the heavy quark approach does not seem to be suit- 
able for the strange sector. The binding energies simply turned out too low for reasonable 
predictions of the mass differences between the heavy baryons and the nucleon. 

On the other hand an interesting path to pursue would be the extension of the light 
sector to flavor SU(3). This would make possible the description of baryons like S c or 
H£. This would in particular be interesting for the issue of flavor symmetry breaking ||. 
Unfortunately the vector meson model for three flavors requires the introduction of additional 



induced components for the strange degrees of freedom [30 
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Appendix A: Bound State Lagrangian 



In this appendix we present the Lagrangian for the ansatze ( j2.13] )-( 2.16| ) of the bound 



heavy mesons. These expressions have already been presented in appendix A of ref [16 



Unfortunately some typographical errors have occurred in the formulas reported there. It is 
therefore appropriate to list the corrected expressions. The present notation corresponds to 
eq O. 

Substituting ((2.15 ) and ( 2T6 ) in ( 2.5|) gives for the S-wave channel 



J(S) 



+- 

r 

-2 



drr 2 y&' 2 + 
2 



M 



a 



R l 
2r 2 



$ 2 + M 



*2 



Rn 



i 2 



+2d{ F' 



a 
— < 
2 

1 

r 



2 

Rn 



f 2 



d2 



9' 



+ 2^2 - * 

1 2 



a 



2r 

1 + cosF) ^o*2 



+ 2Md 

a 



F'ty 1 + -sinF^ 2 



UJ 



+ -sinF 

r 



^2% 



4V2cM 
gm v 



2G' 



a \ R a + 2 

-UJ ^!^ 2 H ^0^1 

2 J 1 2 2r 01 



G 



W^o^l + ^1*2 + -ir (G + 2) * 
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1 V ~ 2r ^ U '- R a W 2 + ^ (G + 2) M> $' + — R a m 



gray I r 



a 



9C C 



(A.1) 



Here a prime indicates a derivative with respect to the radial coordinate r. Furthermore the 
abbreviation R a = cosF — l + + G — cosF) has again been used. 

For the P-wave channel one obtains upon substitution of the ansatz ( |2.13| ) and ( |2.14| ) 

a \ 2 2 ' 1 ^ 



j(P) 



drr 2 1 <l> /2 + 



M — e 



1 + 



1 

+ 2 



tf' 2 - -* 2 

r 



$2 _j_ j^f*2 
1 



*i + 2*2 - ^0 



2r 2 



a 



2 I 
" 2 



Ra - / (i 



— d< -sinF 

r 



-*o + (e - -u ) * 2 



a 



i 2 



+- 



F' 



a 
— < 
2 



e - -uo ) 



1 - cosF) 



2Md 



_,, _ sinF 
F'^i ^ 2 



2\[2cM 

gm v 
A\f2c f 1 



2C C 
2a/*„*i - — W 2 + y~2 (G + 2) ^ 



, . e - ^ ) [G (G + 2) *i - rG'^ 2 ] $ - - 
gmy [r z V 2 / r 

+4 [G (G + 2) $' + G' (2 + F a ) $] V \) • 



' 1 + ^ 



(A.2) 



The typographical errors in ref |L6[ only affect the expressions involving the parameter c. 

Appendix B: Hyperfine Parameters 

In this appendix we give the explicit expressions for the hyperfine splitting parameters 
used in section 4. For convenience we employ additional abbreviations with regard to the 
light meson profiles defined in eqs ( |2.4[) and (|3.3| ) 

Vt = cosF - a (£i - 1 + cosF) , 
V 2 = i-o(& + 6) • 

The explicit expression for the P-wave hyperfine parameter, which enters the mass formula 
for the even parity heavy baryon ( |3.11|) , reads 

2 r°° 



Xp = -j Q drr 2 p[ P H 



(B.l) 
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a 



2a 



e _ u \{y 2 - 2V 1 ) - — {2 + R a )ip 



+ (2Vi + V 2 ) 



a 



4a 

2a 



e - -a; *2 H 



+— (p$ x % + -* 2 + — ^! - — (2 + + 4Md sinF$$ 



d . 
— < smF 



4a 

{2 + R a + VI) ^1^2 ^o^i 



+ F' 



-V 2 ^l + 2a^^ ^2 



| (3£ + £) *o*i + ^ (2 - 26 - 6) ^o*2 + ^'W 2 + ^¥># 



4\/2cM 

4v/2c 
grrip 

+- (2 + # Q ) (26 + 6-2) + 4 [2^' + (2 + i2a) <^'$ - aG'tp®} ^ 



e - ^a; ) ( 4v>*i + ) * + ( V x - ^ 



a 



^ (G + 2) ^! - — V|>2 



(G + 2) 6$' + ( 1 + 2 R <* ) (£l + Q $ - «^V$ 



(B.2) 



It can easily be verified that the terms involving — ^lo^J m M cancel when the heavy 

limit relations for the radial functions ( |2.17| ) are substituted. Taking into account that 
the radial functions which parametrize the heavy meson wave-functions are normalized to 
1/J\e\ « \j\fM (c/. eq ( |2.21| )), it is obvious that the hyperfine parameter \p vanishes in 
the heavy quark limit.. 

The hyperfine parameter for the odd parity baryon (cf. eq ( |3.13|) ) is found to be 



Xs 



2 r°° 

3 Jo 



dr r 2 p^(r) 



a 



2a 



e--iu) (2Vx + V 2 ) + ^rR a <f 



+2 



2 

R n + 2 



$ 2 + (2V1 - V 2 ) 



e u *i 

2 ; 



(B.3) 



2r 



e - |u; ) * 2 



V 2 ^2 + 2-^o + 4—R a (f^ 2 



a 



+4— ^ 



T / t Ra +2 

r^' 2 + ^ 2 ^ — *i 



4a 



- AMd sinF$^ c 



d { F' ( V 2 ^ 2 + — <^ *2 + -sinF 



2a 



(Vl + E a ) ^ 2 + — ^q 

r 



4V2cM r , , 



m v g 



{(£2 " &) *o*i + - (G + 2) 6W2 + -y»'* 1 * 2 + ^^2} 



( -G (26 +6-2) $'vf 2 + % (G + 2) 6<t>*i 



+- [Ra (£ + £) + 2awV] + -5 [2^' - (iW - «G» $] 
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a 



+— [e - -u) (<^i + r(p'y 2 ) $ 



\g (G + 2) + 



-\v 1 + -v 2 



(B.4) 



Again, it can easily be verified that the terms involving (e — ^u?j M vanish when the 
heavy limit relations for the radial functions ( 2.18Q are substituted. With regard to the 
normalization condition ( |2.22| ) the S-wave hyperfine parameter also behaves like xs ~ 1/M 
in the heavy quark limit. 
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